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THEORIES  OF  HIGH  LATITUDE  IONOSPHERE  IRREGULARITIES: 

A  REVIEW* 


1.  INTRODUCTION 

Using  a  variety  of  experimental  techniques,  e.g.,  satellites  f Dyson , 
1969;  Dyson  et  al  ♦ ,  1974;  Sagalyn  et  al.,  1974;  Clark  and  Raitt,  1976;  Phelps 
and  Sagalyn,  1976;  Rodriguez  et  al.,  1981],  rockets  [Olesen  et  al . ,  1976; 
Ogawa  et  al.,  1976;  Kelley  et  al . ,  1980],  scintillations  [Aarons  et  al.,  1969; 
Fremouw  et  al  ♦ ,  1977;  Erukhimov  et  al.,  1981],  and  radar  backscatter  [Weaver, 
1965;  Greenwald , ,  1974;  Hower  et  al.,  1966;  Vickrey  et  al . ,  1980;  Hanuise  et 
al . ,  1981],  it  is  now  known  that  the  high  latitude  ionosphere,  from  the 
auroral  zone  into  the  polar  cap,  is  a  highly  structured  and  nonequilibrium 
medium  containing  Irregularities  (plasma  density  fluctuations  and  structures) 
with  scale  sizes  ranging  from  hundreds  of  kilometers  to  meters.  Aside  from 
being  an  interesting  scientific  phenomenon,  ionospheric  irregularities  are  of 
practical  Interest  to  the  radio-physics  community  since  they  can  disrupt 
transionospheric  radio  wave  communications  channels  (see  recent  review  by 
Davies  [1981]  and  references  therein). 

Several  theories,  *»  g.,  particle  precipitation,  plasma  instabilities  and 
processes,  and  neutral  iluld  dynamics  have  been  proposed  to  account  for  high 
latitude  ionospheric  irregularities.  Recently,  considerable  quantitative 
progress  has  been  made,  especially  in  the  area  of  ionospheric  plasma 
instabilities,  in  identifying  the  physical  processes  that  can  lead  to  high 
latitude  irregularities.  With  the  advent  of  the  EISCAT  incoherent  scatter 
radar  and  new  radar  facilities  to  be  constructed  in  Greenland,  it  seems  timely 
to  review  the  current  state  of  theoretical  research  dealing  with  high  latitude 
ionospheric  irregularities.  As  a  result,  we  present,  in  this  paper,  an 
overview  of  theories  of  naturally  occurring  high  latitude  E  and  F  region 
ionospheric  irregularities.  (For  a  review  of  artificially _  induced 
irregularities,  e.g.,  from  ionospheric  'heating'  experiments,  see  Fejer 
[1979].)  Emphasis  will  be  placed  on  recent  results,  particularly  with  regard 
to  sources  of  high  latitude  irregularities.  Physical  mechanisms  will  be 
stressed  with  extensive  mathematical  analysis  avoided.  Expertise  in  plasma 

physics  has  not  been  assumed.  In — flection — 2G  .we  summarize  the  recent  .and 

"  Or  1 

previous  theoretical  literature  dealing  with  F  region  irregularities^  while. „1  q- 
'Section  (i  we  concentrate  on  E  region  irregularity  phenomena.  . 

.  ■  i. 

•Text  of  invited  talk  presented  at  URSI  Symposium  on  “Radio  Probing  of  the  High-Latitude  Ionosphere 
and  Atmosphere,"  Fairbanks,  Alaska,  August  1982. 
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2.  HIGH  LATITUDE  F  REGION  IRREGULARITIES 


The  high  latitude  latitude  ionosphere  is  profoundly  affected  by  particle 
precipitation,  currents  both  parallel  and  perpendicular  to  the  geomagnetic 
field,  plasma  transport,  and  thermospheric  heating.  In  the  high  latitude  F 
region  ionosphere,  several  irregularity  source  mechanisms  have  been 
proposed:  particle  precipitation,  plasma  processes  and  instabilities,  and 
neutral  fluid  dynamics.  We  now  discuss  each  of  these  sources  separately  and 
point  out  the  irregularity  spatial  scale  size  regimes  in  which  they  are 
believed  to  operate. 

2.1  Structure  from  particle  precipitation 

Particle  precipitation  is  expected  to  play  an  important  role  in 
structuring  the  high  latitude  F  region  ionosphere.  Rees  [1963]  has  shown  that 
low  energy  10^-10®  eV  electrons  deposit  most  of  their  energy  at  auroral  F 
layer  altitudes.  As  a  result,  spatial  and  temporal  structure  in  the  electron 
fluxes  might  be  mirrored  in  the  ambient  ionospheric  plasma.  Evidence  for  this 
hypothesis  was  presented  by  Dyson  and  Wlnnlngham  [1974]  who  showed  good 
correlation  between  structure  in  low  energy  electron  fluxes  and  electron 
density  in  the  polar  cusp.  In  addition  they  indicated  that  the  equatorward 
boundary  of  irregularities  in  the  pre-noon  cusp  is  nearly  collocated  with  the 
boundary  for  soft'  (~  300  ev)  auroral  electron  fluxes.  The  structure  inherent 
in  the  electron  fluxes  presumably  derives  from  processes  occurring  in  the 
plasma  sheet  but  no  quantitative  studies  of  this  structure  have  been  performed 
to  date.  Recently,  Kelley  et  al.  [1982]  gave  further  evidence  that  structured 
low  energy  electron  precipitation  is  a  source  of  large  scale  (X  >  10  km)  high 
latitude  F  region  ionization  irregularities.  They  performed  a  spatial  Fourier 
analysis  of  the  electron  density  irregularities  found  by  Dyson  and  Wlnnlngham 
[1974]  and -found  irregularity  scale  sizes  distributed  from  approximately  75  km 
to  0.75  km.  The  power  spectra  of  these  irregularities  could  be  described  by  a 
power  law  proportional  to  k”**®^. 

Recently,  Vickrey  et  al.  [1980],  using  the  Chatanika  Incoherent  scatter 
radar,  have  studied  large  scale  aproxlmately  magnetic  field  aligned  convecting 
plasma  enhancements  in  the  midnight  sector  auroral  F  region  (see  Fig.  1). 


These  ionization  enhancements,  which  can  have  their  plasma  density  enhanced  by 
factors  of  two  to  ten  over  background  ambient  values,  have  been  observed  in 
regions  of  diffuse  auroral  particle  precipitation  and  associated  field  aligned 
currents.  The  occurrence  of  the  plasma  enhancements  is  apparently  not 
strongly  related  to  magnetic  activity  or  to  E  region  processes.  The  overall 
north-south  dimensions  of  some  of  the  observed  plasma  enhancements  is 
comparable  to  the  outer  scale  sizes  of  the  electron  density  structures 
associated  with  auroral  F  region  particle  precipitation  events  [Kelley  et  al., 
1982].  As  a  result,  the  structured  low  energy  particle  precipitation  may  be 
an  important  source  of  the  large  scale  F  region  structure  associated  with 
these  plasma  enhancements.  Large  scale  east-west  structure  associated  with 
the  plasma  enhancements  has  also  been  demonstrated  [Tsunoda  and  Vickrey, 
1982].  In  addition,  spatial  spectral  analysis  of  these  enhancements  [Kelley 
et  al.,  1982]  have  indicated  structure  from  approximately  50  km  down  to 
several  kilometers.  The  presence  of  smaller  kilometer  and  hundreds  of  meter 
sized  plasma  density  irregularities  collocated  with  these  enhancements  has 
been  verified  using  satellite  scintillation  studies  [Fremouw  et  al . ,  1977; 
Rlno  et  al.,  1978;  Vickrey  et  al.,  1980].  The  scintillation  measurements  have 
indicated  that  the  electron  density  irregularities  are  structured  like  L  shell 
aligned  sheets  for  equatorward  convection  of  the  plasma  enhancements. 

Similar  observations  [Weber  and  Buchau,  1981]  of  correlations  between 
structured  strongly  enhanced  F  region  plasma  and  low  energy  electron 
precipitation  have  also  been  reported  for  the  polar  cap  ionosphere. 

2.2  Structure  from  plasma  processes  and  instabilities 

It  is  logical  to  consider  plasma  processes  in  discussing  high  latitude 
ionospheric  irregularities  since  several  sources  of  free  energy  are  available 
to  drive  various  plasma  instabilities.  Examples  of  these  sources  include 
density  gradients,  velocity  shears,  and  currents  both  parallel  and 
perpendicular  to  the  geomagnetic  field.  Both  plasma  macroinstabilities,  which 
are  fluidlike  and  operate  on  scale  sizes  A  »  p  ,  p^,  with  PgtPj)  being  the 
electron  (ion)  gyroradius,  and  plasma  microinstabilities  (A  <  pg,  p , )  have 
been  Invoked  to  account  for  high  latitude  ionospheric  irregularities. 
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V  .  J  -  0  (5) 

Here  n^Ca  »  i  or  e)  is  the  species  density  and  ie  the  total  electric 
field.  Since  we  will  be  interested  in  low  frequencies  and  long  wavelenghts, 
we  have  ignored  inertial  terms  in  the  electron  and  ion  momentum  equations  (3) 
and  (4).  Equation  (5)  results  from  the  assumption  of  quasi-neutral 
fluctuations  n  “  n.  ■  n.  In  addition  vrt  and  V-  refer  to  the  electron  and  ion 
velocities  along  the  magnetic  field  giving  rise  to  a  cold  magnetic  field 
aligned  current.  The  symbol  denotes  the  ion-neutral  collision 

frequency,  v  the  electron-ion  collision  frequency,  c  the  speed  of 
light,  T  *  T.  =  T  the  species  temperature,  ce  the  ion  acoustic  speed, 
and  )  the  ion  (electron)  gyrof requency .  We  have  neglected  vgn  compared 

with  v  ^  and  taken  «  1  for  a  ■  i,  e  (F  region  approximation). 
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The  following  geometry  is  used:  the  y  axis  is  in  the  north-south 
direction,  the  x  axis  points  west,  and  the  z  axis  is  downward  along  the 
magnetic  field. 


Any  two  equations  (1),(2),  and  (5)  provide  a  complete  description  of  the 

problem.  We  will  use  the  ion  continuity  equation  (1)  and  (5).  After 

separating  the  total  electric  field  into  an  ambient  and  fluctuating 

part  E.  •  E  —  V .  and  transforming  to  a  frame  drifting  with  velocity  V-  -  - 
—A.  — O  1  — ° 

(c/B)  [£  x  we  can  linearize  (1)  and  (5)  by 

separating  n  «  n  (y)  +  6n  with  6n,  6\J»  «  exp  [i(k  x  +  It  z  -  ut)J, 

_  O  x  z 

u)  •  (i)  +  1y,  L  1  5  (1/n  )(3n  /3y).  This  gives  a  growth  rate 
r  o  o 


v  ,  ,  v.  cE 

*  (®2  +  t r1  '  Dikx2 '  D«kz2 

i  e 


where  0  =  kz/kx,  Vj  -  _z(vQ  -  VQ),  kj2  -  kx2,  k±2  »  kz2,  D±  -  (vei/ne%>cs2 

and  D,  -  (cs2/vin){l+  Uvin/ni)2/<<veivin/Qeni>  +  (k22/ki2))]} *  (•  denotes 

along  z). 


The  expression  for  the  growth  rate  y  in  (6)  can  be  maximized  as  a 

function  of  6  >  k  /k  ,  a  measure  of  field  alignment,  using  Sy/ddle-e  “  0 
Z  X  Dl 

giving 


-4 

Using  typical  diffuse  auroral  F  region  parameters  v.  /ft.  =  10  , 

A  *  e  _  o 

v  /ft  *10  ,  E  3  10  mV/m,  j  ■  n  eV.  “  1  yA/m2,  B  «  0.5  G,  n  *10  cm  , 

this  gives  | 8  |  =  10  ,  i.e.,  approximate  field  alignment.  Inserting  these 

m  8 
parameters  into  (6)  with  L  3  20  km,  *  0.2  m2/s,  and  »  10°  m2/s, 

we  find  that  the  fastest  growing  linear  modes  have  growth 

times  Y  *  *  102  s. 
max 

Although  there  is  only  one  growth  rate  y  as  given  by  (6)  it  is  convenient 

to  discuss  y  in  two  limits.  If  kz  «  0,  then  the  growth  rate 

Y  3  cE  /BL  (E  x  gradient  drift  instability)  while  if  E  ■ 

°  -1  /2 

0,  y  3  (V./2L)(1  +  v,  ft  /ft.v  )  '  and  results  from  the  current 

max  '  d  in  e  i  e' 

convective  instability.  We  now  discuss  these  two  instabilitiles . 
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It  is  well  know  that  convecting  ionospheric  plasma  clouds  are  unstable, 
under  certain  conditions,  and  can  produce  plasma  Irregularities  through  the  E 
x  gradient  drift  instability.  The  gradient  drift  instability  [Simon, 
1963;  Linson  and  Workman,  1970]  is  a  convective  instability  with  its  nonlinear 
evolution  [Zabusky  et  al.,  1973;  Scannapieco  et  al.,  1976]  resembling  the 
classical  Rayleigh-Taylor  instability  [ Chandrasekhar ,  1960]  which  arises  when 
a  heavy  fluid  is  supported  by  a  lighter  fluid.  Ihe  basic  JJxB^  instability 
mechanism  can  be  understood  by  noting  Figure  2.  Here  the  horizontal  line 
represents  an  unperturbed  contour  of  electron  density.  In  addition,  the 
ambient  background  electric  field  is  in  the  y-direction,  the  ambient 
magnetic  in  the  z-direction,  and  the  background  density  gradient  points  in 
the  x-direction.  Let  the  density  be  perturbed  by  a  small  amplitude  sinusoidal 
perturbation  with  wavenumber  Jc  parallel  to  .  In  the  F  region,  the  ions 
drift  to  the  right  in  the  Pedersen  direction  relative  to  the  electrons.  This 
gives  rise  to  space  charges  (+  and  -)  which  in  turn  cause  small  scale  electric 
fields  _E_'  alternating  in  direction  as  shown.  The  corresponding  E*xB0  drifts 
will  then  convect  depleted  regions  upward  (toward  increasing  density)  and 
enhanced  regions  downward  (toward  decreasing  density)  with  the  result  that 
they  both  appear  to  grow  relative  to  the  background  density  gradient  -  an 
unstable  situation.  In  the  previous  configuration  the  E^  x  B  gradient  drift 
instability  can  also  arise  with  a  neutral  wind  _U  blowing  in  the  -x  direction 
with  no  E0 . 

Using  both  analytical  and  numerical  simulation  techniques,  Reskinen  and 
Ossakow  [1982a]  have  studied  the  linear  stability  and  nonlinear  evolution  of 
large  scale  convecting  plasma  enhancements  in  the  auroral  F-region 
ionosphere.  Their  results  indicate  that  convecting  diffuse  auroral  plasma 
enhancements  can  be  driven  unstable  through  the  ExB  gradient  drift 
instability.  This  destabilization  can  both  directly  and  indirectly  generate 
plasma  density  and  electric  field  irregularities  and  fluctuations  with  scale 
sizes  of  several  hundred  kilometers  to  tens  of  meters.  These  irregularities 
take  the  form  of  anisotropic  striation-like  structures  (elongated  in  the 
north-south  direction  for  equatorward  convection)  which  can  form  on  the  order 
of  half  an  hour  under  typical  auroral  conditions.  Fig.  3  illustrates  the 
model  of  the  plasma  enhancement  used  in  the  numerical  simulations.  In  a  plane 
nearly  perpendicular  to  the  geomagnetic  field  Fig.  4  displays  the  linear  and 
nonlinear  evolution  of  the  ExB^  Instability  in  a  convecting  auroral  F  region 
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plasma  enhancement.  Parameters  were  used  that  agree  with  the  observations  of 

Vickrey  et  al.  [1980J.  Fig.  4a  shows  the  initial  configuration  of  the 

Ionization  enhancement  and  includes  a  small  12  density  perturbation.  Fig.  4b 

illustrates  the  linear  regime  at  t  *  550  sec  and  shows  unstable  growth  on  the 

poleward  side  of  the  equatorward  convecting  plasma  enhancement.  Fig.  4c  gives 

the  structure  of  the  plasma  enhancement  at  t  =  1000  sec.  Finally,  Fig.  4d 

displays  the  plasma  enhancement  at  t  -  1600  sec  in  the  we  11 -developed 

nonlinear  regime  where  steepened  and  elongated  striation-like  structures  are 

evident.  These  fingers  or  striations  will  be  oriented  in  a  direction 

dependent  upon  the  local  electric  field  magnitude  and  direction.  These 

theoretical  results  are  not  inconsistent  with  recent  experimental  observations 

[Tsunoda  and  Vickrey,  1982],  using  the  Chatanika  radar,  which  Indicate  large 

neutral  wind  velocities  (_E_  x  gradient  drift  instability)  and  fingerlike 

structures  collocated  with  large  scale  convecting  plasma  enhancements  (see 

Fig.  5).  Vickrey  and  Kelley  (1982]  have  studied  the  role  of  classical 

diffusion  and  a  conducting  E  layer  in  removing  these  irregularities  once  they 

are  produced.  In  addition,  Keskinen  and  Ossakow  [1982a]  found  that  the  larger 

scale  size  irregularities  (fingers)  can  cascade  to  smaller  scale  size 

structures  through  nonlinear  mode  coupling  and  two-step  processes.  Examples 

of  the  spatial  power  spectra  of  these  irregularities  in  the  north-south  P(ky) 

and  east-west  P(kx)  direction  taken  in  the  nonlinear  regime  of  the  simulations 

is  shown  in  Fig.  6.  These  power  spectra  can  be  well  represented  by  power  laws 

-n„ 

of  the  form  P(kx)  «  kx  .  with  n*  -  2-2*5  for  2x/kx  between  100  km  and  3  km 
while  in  the  north-south  direction  P(ky)  *  ky  ^  with  ny  =  2  for  2x/k  between 
256  km  and  3  km.  This  process  of  finger  formation,  elongation,  and  steepening 
is  almost  self-similar  in  character  with  similar  morphologies  and  power 
spectra  for  scale  sizes  X  between  1  km  and  100  m  [Keskinen  and  Ossakow, 
1982b].  Some  observational  evidence  [Vickrey  et  al.,  1980]  indicates  that 
these  plasma  enhancements  are  probably  subjected  to  ambient  auroral  convection 
patterns.  As  a  result,  these  enhancements  could  be  a  major  source  of  F  region 
ionospheric  irregularities  throughout  the  auroral  zone  and  polar  cap. 

The  _ExB^  instability  in  large  scale  plasma  enhancements  results  from  a 
coupling  between  the  convective  electric  field  and  a  density  gradient 
perpendicular  to  the  magnetic  field.  However,  the  coupling  of  density 
gradients  and  magnetic  field  aligned  currents  can  also  lead  to  plasma 
instability  through  the  current  convective  instability  [ Lehnert ,  1958; 


Kadomtsev  and  Nedospasov,  I960].  Ossakow  and  Chaturvedi  [1979]  showed  that 


plasma  enhancements  can  also  be  linearly  unstable  to  the  current  convective 

Instability  in  regions  where  the  _ExJ^  gradient  drift  instability  is  stable, 

i.e.,  in  regions  where  the  local  convection  velocity  ?<nd  density  gradient  are 

in  opposite  directions.  The  basic  physical  mechanism  responsive  for  the 

current  convective  instability  is  as  follows  (see  Fig.  7):  Let  the  magnetic 

field  and  current  jQ|  be  in  the  z-direction,  an  ambient  electric  field  Js^ 

with  perpendicular  component  in  the  x-direction,  and  a  density  gradient  in  the 

y-direction.  Consider  a  density  perturbation  with  wavevector  as  shown.  The 

projection  on  Jc  of  the  ion  Pedersen  drift  caused  by  E  ^  results  in  a  drift 

that  is  stable  to  the  j&B  instability.  However,  the  assumed  direction  of 

implies  a  relative  drift  between  ions  and  electrons  (in  the  electron  rest 

frame)  is  anti-parallel  to  Bn .  This  motion  projected  on  Jc  gives  space  charges 

and  subsequent  electric  fields  E” .  If  the  particle  motion  projected  on  k  is 

dominated  by  J  then  the  total  perturbation  space  charge  electric  fields  will 

be  denoted  by  E' .  The  corresponding  E'  x  ^  drifts  will  then  convect  enhanced 

(depleted)  regions  out  of  (into)  the  figure  which  will  appear  to  grow  relative 

to  the  background  density  in  direct  analogy  to  the  ExB^  gradient  drift 

instability  as  outlined  previously.  The  ratio  of  the  linear  growth  rate  of 

the  £  x  B  gradient  drift  instability  versus  the  current  convective  instability 

can  be  written  approximately  (2V  /V,)(l  +  ft  v./ft. v  )^^  »  2V  /V,  *  2V  ne/j  , 

o  a  el  l  e  o  a  o  ol 

since  BgVj/  ftjVe  is  of  order  unity  in  the  auroral  F  region.  Here  VQ  -  cEQ/B 
is  the  convection  speed  across  the  magnetic  field.  As  a  result,  for  strong 
perpendicular  electric  fields  and/or  weak  currents  and  ambient  densities  the 
ExB  gradient  drift  instability  will  dominate  and  vice  versa.  The  nonlinear 
evolution  of  the  current  convective  instability  in  the  auroral  F  region 
ionosphere  was  investigated  by  Kesklnen  et  al .  [1980].  The  evolution  and 
morphology  of  the  plasma  density  irregularities  generated  by  the  current 
convective  instability  were  similar  to  that  of  the  ExE^  instability. 
Basically,  the  current  convective  instability  saturates  (stabilizes) 
nonlinearly  by  feeding  energy  from  the  linearly  unstable  waves  to  the  linearly 
damped  harmonics  [Chaturvedi  and  Ossakow,  1979].  Recently,  the  current 
convective  instability  in  the  auroral  ionosphere  has  been  extended  to  include 
other  effects,  e.g.,  magnetic  shear  [Huba  and  Ossakow,  1980]  and  ion  inertial 
and  collisional  effects  [Chaturvedi  and  Ossakow,  1981]. 
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Another  source  of  free  energy  In  the  high  latitude  F  region  ionosphere  is 
the  magnetic  field  aligned  current  system.  Drummond  and  Rosenbluth  [1962] 
first  outlined  the  theory  of  field  aligned  current  driven  ion  cyclotron 
instabilities  for  a  collisionless  plasma  (for  k^  >  1).  Kindel  and  Kennel 
[1971]  applied  this  theory  to  the  auroral  zone  on  the  topside  ionosphere  where 
the  plasma  is  only  weakly  collisional.  Chaturvedi  [1976],  by  invoking  the 
collisional  ion  cyclotron  instability,  has  shown  that  strong  field  aligned 
currents  in  the  auroral  F  region  ionosphere  can  lead  to  irregularities  with 
scale  sizes  of  hundreds  of  meters  perpendicular  to  the  magnetic 
field  (kj^Pj  <  1).  The  basic  mechanism  is  that  a  field  aligned  current  can 
excite,  owing  to  dissipative  effects  due  to  electron-neutral  collisions,  a 
growing  ion  cyclotron  wave  propagating  nearly  transverse  to  the  geomagnetic 
field.  Current  velocities  of  several  kilometers  per  second  are  needed  to 
excite  this  instability. 

Hudson  and  Kelley  [1976]  have  demonstrated  that  a  temperature  gradient 
driven  drift  wave  instability,  which  arises  from  colllnear  density  and 
temperature  gradients  perpendicular  to  the  magnetic  field,  might  explain 
density  irregularities  with  scale  sizes  of  hundreds  of  meters  observed  at  the 
equatorward  edge  of  the  ionospheric  plasma  trough. 

Kelley  and  Klntner  [1978]  have  argued  that  highly  structured  electric 
fields,  presumably  of  magnetospheric  origin  commonly  seen  in  the  dayslde 
winter  high  latitude  ionosphere,  might  cause  density  irregularities  with 
similar  scale  sizes  through  the  mixing  of  flux  tubes  that  have  varying  plasma 
density. 


2.2.2.  Plasma  microinstabilities 

Small  scale  irregularities  with  sizes  on  the  order  of  and  smaller  than 
the  ion  gyroradius  (~  10  m  in  the  auroral  F  region)  can  be  generated  very 
efficiently  through  various  plasma  instability  processes.  These  small  scale 
irregularities  could  probably  be  detected  by  most  high  latitude  backscatter 
radars . 

Plasma  density  gradients,  both  sharp  and  smooth,  are  a  ubiquitous  feature 
of  the  high  latitude  F  region  ionosphere.  These  density  gradients  can  drive  a 
variety  of  drift  wave  type  plasma  Instabilities.  If  is  the  ion  Larmor 
radius  and  L  the  density  gradient  scale  length,  then,  in  an  approximate  sense, 
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the  following  modes  are  expected  to  dominate  in  the  following  density  gradient 

1/2 

regimes  [Mlkhallovskil ,  1974]  p^/L  <  (me/m^)  ,  universal  drift  mode; 

(mg/mj)1^2  <  Pj/L  <  (me/mi)1^,  drift  cyclotron  mode;  (nig/m^)1^  <  pj/L, 

lower-hybrid-drift  mode.  The  growth  of  these  drift  modes  tends  to  maximize 
for  wavelengths  X  near  the  electron  or  ion  Larmor  radius.  Basu  [1978]  has 
noted,  using  Ogo  6  satellite  observations,  a  high  correlation  between  small 
irregularity  structures  and  high  latitude  sub-trough  plasma  density  gradients 
and  invoked  drift  waves  to  explain  these  fluctuations.  Other  sources  of  sharp 
density  gradients  in  the  high  latitude  F  region  can  be  found  near  auroral  arcs 
and  near  large  scale  convecting  plasma  enhancements  [Vickrey  et  al . ,  1980]. 

Another  source  of  plasma  free  energy  that  could  directly  excite  small  and 
large  scale  irregularities  in  the  high  latitude  F  region  ionosphere  is 
velocity  sheared  plasma  flows,  e.g.,  near  auroral  arcs  [Kelley  and  Carlson, 
1971].  For  unstable  velocity  sheared  plasma  flows  perpendicular  to  a  magnetic 
field,  the  growth  is  nonlocal  and  maximizes  for  irregularity  wavenumbers  k 
such  that  kLy  <  1  [Mlkhailovskli ,  1974]  where  Ly  is  the  velocity  shear 

gradient  scale  length.  A  local  instability  (kl^  >  1)  can  operate  for  velocity 
sheared  plasma  flow  parallel  to  a  magnetic  field  [Mikhallovskll ,  1974]. 

Recently,  Keskinen  and  Huba  [1982],  using  kinetic  theory,  showed  that  velocity 
sheared  electron  flow  parallel  to  the  geomagnetic  field  near  strong  discrete 
auroral  arcs  might  also  be  unstable  and  lead  to  centimeter  sized  small  scale 
irregularities  with  scale  sirfes  on  the  order  of  the  electron  gyroradius. 

Velocity  space'  instabilities,  as  opposed  to  the  previously  discussed 
configuration  space  instabilities,  have  also  been  invoked  to  explain  high 
latitude  F  region  small  scale  irregularities.  Ott  and  Farley  [1975]  have 

found  that  the  action  of  ion-neutral  charge  exchange  collisions  can  lead  to 
anisotropic  F  region  ion  velocity  distributions  under  the  influence  of  large 
ExB^  convection  velocities.  They  showed  that  such  distributions  are  unstable 
to  the  Post-Rosenbluth  instability  [Rosenbluth  and  Post,  1965;  Post  and 
Rosenbluth,  1965]  which  is  a  short  wavelength  (k  »  fl^/v^),  high  frequency 
(u)  »  Jip  instability  in  which  is  nearly  perpendicular  to  Here 

k,  a),  v^  are  the  wavenumber,  frequency,  ion  cyclotron  frequency,  and  JJxB^ 
drift  speed,  respectively.  They  showed  that  this  instability  would  be  excited 
at  short  wavelengths  (10-20  cm)  and  argued  that  it  might  lead  to  density 
fluctuations  of  a  few  percent. 
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2.3  Structure  from  neutral  fluid  dynamics 


Since  the  high  latitude  ionosphere  is  a  coupled  medium  consisting  of  both 
plasma  and  neutral  constituents,  various  plasma  processes  and  irregularities 
could  be  transferred  to  the  neutral  gas  via  colllsional  effects  and  vice 
versa.  Using  the  Chatanika  radar,  several  studies  [ Baron,  1972;  see  recent 
review  by  Hunsucker ,  1982  and  references  therein]  have  shown  that  the  high 
latitude  ionosphere  often  exhibits  quasi-periodic  fluctuations  in  electron 
density  and  temperature  with  these  fluctuations  having  large  equatorward 
velocities.  It  is  generally  accepted  [Hunsucker,  1982]  that  these  traveling 
ionospheric  disturbances  [TID's)  are  ionospheric  manifestations  of  atmospheric 
gravity  waves  [Hines,  I960].  The  TID's  are  a  very  large  scale  phenomena  with 
wavelengths  of  several  hundred  to  thousands  of  kilometers.  The  most  likely 
sources  [Hunsucker,  1982]  of  the  gravity  waves/TID's  are:  (1)  Joule  heating 
and  Lorentz  forces  associated  with  the  auroral  electro jet  and  (2)  intense 
particle  precipitation  events.  These  large  scale  TID's  could  act  to  seed 
other  macroscopic  plasma  Instabilities.  Much  work  remains  to  be  done  in  this 
area. 
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3.  HIGH  LATITUDE  E  REGION  IRREGULARITIES 


In  Che  high  latitude  E  region,  experimental  observations  using  both  radar 
backscatter  fHarang  and  Stoff regen,  1938;  Bowles ,  1955;  McNamara ,  1955; 
Balsley  and  Ecklund,  1972  see  reviews  by  Hultqvist  and  Egeland  (1964),  Lange- 


Hesse  (1967),  Unwin  and  Baggaley,  (1972)  Greenwald  (1979),  Fejer  and  Kelley 
(1980)]  and  rocket  techniques  [Kelley  and  Mo2er,  1973;  Holtet ,  1973;  Olesen  et 
al . ,  1976;  Ogawa  et  al . ,  1976]  have  shown  the  existence  of  plasma  density 
irregularities  with  scale  sizes  of  meters  to  hundreds  of  meters.  These 
irregularities  have  been  shown  to  be  highly  correlated  with  the  auroral 
electro jet  [Greenwald  et  al.,  19  75b;  Tsunoda  et  al  ♦ ,  1976a].  Various  plasma 
instabilities  and  processes  have  been  successful  in  accounting  for  many  of  the 
features  of  these  irregularities. 

Fluid  equations  that  can  describe  the  basic  dynamics  of  the  auroral  E 
region  electro jet  plasma  can  be  written  [Rogister  and  D'Angelo,  1970;  Sudan  et 
al.  1973] 
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where  E  ■  -  V<|»  and  n  ,  v  ,  v  ,  m  ,  T  are  (a  =  e,i)  the  density,  velocity, 
—  a  —a  a  a  a 

collision  frequency  with  neutrals,  mass,  and  temperature,  respectively.  In 

addition,  ip  is  the  electrostatic  potential,  E  the  electric  field,  c  the  speed 

of  light  and  ji  the  magnetic  field.  Equations  (8)  -  (12)  are  valid  under  the 

assumptions  of  electrostatic  waves,  charge  neutrality  ng  =  ,  isothermal 

electrons  and  ions,  and  neglect  of  electron  inertia.  Taking  B  « 

B.X,  E  -  e  x  -  V60(  linearizing  eq.  (7)  -  (12)  with 
o—  —  o— 

6n^,  6<V  *  exp[(J<  •  3^  ~  wt)]  and  assuming  a  relative  drift  between 
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electrons  and  ions  across  the  magnetic  field,  we  find  the  frequency  and  growth 
rate  of  waves  propagating  perpendicular  to  the  magnetic  field: 

\  -  k  •  Vj/d  +  *)  (13) 

Yk  (♦/l-N»)[(fle/ve)(wk/kL)  +  (wk2  -  k2cg2)/v1J  (14) 

where  \]»  ■  v^v  ^  *  In  /Vn  I  with  Vn  the  electron  density  gradient  along 

E^,  and  Cg  =  [(K(Te  +  T^)/m^]^^  the  ion  acoustic  velocity.  Although  (13)  and 
(14)  describe  only  one  instability,  it  is  convenient  to  discuss  two  well  known 
limits.  If  L  ♦  «°  (no  density  gradient)  then  (13)  and  (14)  reduce  to  the 
Farley-Buneman  instability  [Farley ,  1963;  Buneman ,  1963]  with  instability 

threshold  condition  >  cg.  If  «  cg  and  ( /J2k ) (  1/kL)  >  1  then 
the  u^2  term  in  (14)  can  be  neglected  recovering  the  _E  x  gradient  drift 
(cross  field)  instability  [Simon,  1963]  with  approximate  instability 
criterion  (8e/v  )(V./L)  >  (k2cg2/v^)  when  Jt  is  along  the  y  direction. 

Radar  backscatter  studies  have  provided  much  information  regarding  small 
scale  irregularities  (with  meter-sized  wavelengths)  in  the  high  latitude  E 
region.  The  E  region  radar  echoes,  which  are  highly  aspect  sensitive,  can 
only  be  observed  at  large  angles  to  the  electrojet  current  due  to  geometric 
considerations.  Basically,  the  echoes  can  be  divided  into  two  types  depending 
upon  their  duration  and  latitudinal  and  longitudinal  extent.  Additional 
characteristics  such  as  location,  threshold,  doppler  feature,  etc.  have  been 
summarized  previously  [Hultqvist  and  Egeland,  1964;  Lange -Hesse,  1967;  Unwin 
and  Baggeley,  1972;  Greenwald ,  1979;  Fejer  and  Kelley,  1980].  The  Farley- 
Buneman  (two-stream)  instability,  which  is  driven  by  the  electrojet  current, 
has  been  invoked  to  explain  the  small  scale  irregularities  in  the  auroral  and 
polar  cap  E  region  ionosphere  [Olesen,  1972;  Moore roft ,  1972;  Wang  and 

Tsunoda ,  1975;  Olesen  et  al . ,  1975,  1976;  PrimdaHl  et  al.,  1974].  The  linear 
theory  of  the  Farley-Buneman  instability  seems  to  explain  several  features  of 
small  scale  irregularities  in  the  high  latitude  E  region  especially  in  the 
polar  cap,  e.g.,  peak  backscatter  and  electric  field  fluctuations  in  the 
direction  of  the  current  [Olesen  et  al.,  1976;  Tsunoda  et  al.,  1976a;  Bahnsen 
et  al .  1978].  However,  several  irregularity  characteristics  cannot  be 

explained  using  linear  theory,  e.g.,  waves  propagating  with  small  and  large 
aspect  angles  perpendicular  to  the  electro jet  current,  wave  phase  velocity  in 
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excess  of  the  ion  sound  speed,  saturated  wave  amplitudes  and  spectra,  and 
electron  heating.  Several  nonlinear  saturation  theories  of  the  Farley-Buneman 
instability  have  been  proposed,  e.g>,  quasilinear  effects  [Sato,  1972; 
Rogister ,  1971],  resonance  broadening  [welnstock  and  Sleeper,  19/2],  mode 

coupling  [Rogister  and  Jamin,  1975],  and  stabilization  by  external  low 
frequency  E  x  turbulence  [ Kesklnen,  1981].  These  nonlinear  theories  predict 
different  saturated  wave  amplitudes  and  steady  state  spectra.  To  explain  the 
large  aspect  angle  waves  Hofstee  and  Forsyth  [1972]  and  Moorcroft  [1972]  have 
suggested  an  ambient  magnetic  field  distortion  caused  by  large  electrojet 
currents  while  Volosevich  and  Llperovsky  [1975]  argued  for  ion  acoustic  waves 
generated  by  E  region  field  aligned  currents.  Recently,  St. -Maurice  et  al . 
[1981]  and  Schlegel  and  St. -Maurice  [1981]  have  shown  that  anomalous  electron 
temperatures  in  the  high  latitude  E  region  can  be  quantitatively  explained  in 
terms  of  plasma  heating  caused  by  unstable  Farley-Buneman  waves  (see  Fig. 
18).  Good  agreement  is  made  between  theory  and  observations.  Finally, 
D'Angelo  [1973]  has  proposed  the  ion  cyclotron  instability,  as  generated  by 
field  aligned  currents  in  the  high  latitude  E  region,  to  be  another  source  of 
both  large  and  small  scale  irregularities.  D'Angelo  notes  that  the  ion 
cyclotron  wave  excitation  would  probably  take  place  at  altitudes  higher  than 
those  of  the  auroral  elctrojet  current  and,  subsequently,  propagate  downward 
at  an  a.pect  angle  on  the  order  of  10%. 

Larger  scale  .  irregularities  in  the  high  latitude  E  region  with 
wavelengths  up  to  several  hundreds  of  meters  have  been  identified  principally 
through  in  situ  rocket  probes  [Kelley  and  Mozer,  1973;  Holtet ,  1973;  Ogawa  et 
al . ,  1976].  Vertical  density  gradients  are  often  observed  in  conjunction  with 
auroral  electro jet  currents.  As  a  result,  the  JJxB^  gradient  drift  instability 
has  been  invoked  [Greenwald ,  1974,  1975a]  to  explain  the  large  scale 

irregularities.  This  instability  is  different  from  the  ExB  instability  used 
to  describe  high  latitude  F  region  irregularities  (see  Section  2).  In  the  F 
region  case,  the  basic  current  flow  is  in  the  Pedersen  direction  while  the 
Hall  current  dominates  in  the  E  region.  However,  the  physical  mechanisms 
describing  the  _ExJJ^  instability  in  both  regions  are  essentially  identical. 
Again,  the  linear  theory  of  the  J!xB  instability  cannot  account  for 
irregularities  propagating  perpendicular  to  the  electrojet  current  and 
saturated  wave  amplitudes  and  spectra.  Several  nonlinear  theories  of  the  E- 
region  gradient  drift  instability,  e.g.,  quasilinear  effects  [Sato  and  Ogawa, 


1976],  resonance  broadening  [Welnstock  and  Williams  1971]  and  strongly 
turbulent  nonlinear  wave-wave  Interactions  [Sudan  and  Kesklnen,  1977]  have 
been  proposed.  Greenwald  [1974]  has  attempted  to  account  for  small  scale 
oblique  waves  by  appealing  to  a  two-step  multilinear  process  [Sudan  et  al . , 
1973]  In  which  secondary  vertical  short  wavelength  waves  can  grow  at  the 
expense  of  the  perturbation  electric  fields  and  density  gradients  associated 
with  the  linearly  unstable  primary  horizontal  longer  wavelength  gradient  drift 
fluctuations.  Greenwald's  analysis  shows  that  secondary  gradient-drift  waves 
can  be  excited  while  there  is  difficulty  in  generating  secondary  Farley- 
Buneman  waves.  Radar  echoes  resembling  secondary  Farley-Buneman  modes  have 


4.  SUMMARY 


We  have  attempted  to  give  an  overview  of  theoretical  interpretations  of 
naturally  occurring  high  latitude  E  and  F  region  ionospheric  irregularities. 
Much  quantitative  progress  has  been  made  in  recent  years  in  identifying  the 
source  mechanisms  of  high  latitude  irregularities.  Plasma  instabilities  have 
been  shown  to  play  an  Important  role  in  the  generation  and  evolution  of  these 
Irregularities . 

In  the  auroral  F  region,  there  is  some  evidence  for  large  scale 
(\  >  10  km)  structure  and  irregularities  being  produced  by  particle 
precipitation  events  and  subsequent  formation  of  large  scale  size  plasma 
enhancements.  Some  experimental  observations  indicate  that  these  plasma 
enhancements  appear  to  follow  ambient  auroral  convection  patterns.  Ihese 

convecting  plasma  enhancements  have  been  shown  to  be  unstable  to  plasma 
macroinstabilities,  i.e.,  the  ExB  and  current-convective  instability,  which 
can  generate  irregularities  and  structure  with  scale  sizes  on  the  order  of  and 
smaller  than  the  approximate  size  of  the  plasma  enhancements.  Nonlinear  wave- 
wave  interactions  and  associated  wavenumber  cascades  can  lead  to  further 
structuring  until  diffusive-like  processes  dominate  and  irregularity  formation 
is  quenched.  Since  these  plasma  macroinstabilities  become  highly  nonlinear 
further  progress  in  this  area  can  only  be  made  through  numerical  modeling  and 
simulations  of  the  fundamental  plasma  equations.  Moreover,  since  the 
convective  motions  can  be  spatially  dependent  (e.g.,  ^  can  be  a  function  of 
space)  and  so  introduce  velocity  shear,  it  is  important  to  assess  the 
influence  of  this  velocity  shear  on  the  plasma  macroinstabilities  (see  for 
example  Guzdar  et  al.,  1982;  Huba  et  al.,  1982).  This  would  also  include 

being  able  to  model  this  shear  in  the  nonlinear  numerical  simulations.  In 
addition,  the  role  of  neutral  fluid  disturbances,  i.e.,  gravity  waves,  in 
producing  large  scale  irregularities  should  be  further  Investigated.  These 
waves  could  act  to  "seed"  large  scale  plasma  macroinstabilities.  The 

generation  of  smaller  scale  irregularities  with  the  wavelengths  of  meters  to 
tens  of  meters  in  the  auroral  F  region  has  been  given  considerably  less 
theoretical  attention.  Free  energy  sources  for  these  smaller  scale  structures 
include  sharp  density  gradients  (leading  to  drift  waves),  velocity  sheared 
plasma  flows  both  parallel  and  perpendicular  to  the  magnetic  field,  and 
anisotropic  plasma  velocity  distributions.  An  outstanding  question  is  the 
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saturated  amplitudes  and  spectra  of  the  small  scale  irregularities  and  the 
role  they  play  in  influencing  larger  scale  irregularities  and  structures. 

In  the  high  latitude  E  region,  the  strong  electro jet  currents  across  the 
geomagnetic  field  have  been  shown  to  drive  plasma  instabilities,  i.e.,  the 
Farley-Buneman  and  JjxB  instabilities,  which  can  lead  to  both  small  and  large 
scale  irregularities.  The  linear  theories  of  these  instabilities  can  explain 
several  features  of  the  irregularities.  Results  from  nonlinear  studies  e.g., 
saturated  amplitudes,  of  these  instabilities  are  conflicting  in  several 
respects.  Further  analytical  and  numerical  studies  of  these  instabilities 
especially  the  Farley-Buneman  (two-stream)  instability  are  needed. 
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Fir.  3  —  Model  of  plasma  enhancement  used  in  numerical  simulations  [Keskincn  and 
Ossakow,  1982],  The  quantity  nQ(x,y)  represents  the  initial  plasma  enhancement 
density  profile  while  N0  is  the  average  background  auroral  F  region  plasma  density. 
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Fig.  4  —  Real  apace  plasma  isodensity  contour  plots  of  [n0(x,y)  +  6n(x,y)]/N0  at  (a)  t  ■  0  sec, 
(b)  t  -  660  sec,  (c)  t  ■  1000  sec,  (d)  t  *  1600  sec.  The  y-axis  is  compressed  by  a  factor  of  2.68. 
The  distance  between  tick  marks  in  the  x  direction  (y  direction)  is  6  km  (12.8  km).  Eight  contours 
are  plotted  in  equal  increments  of  1.26  beginning  at  1.26.  The  observer  is  looking  upward  along  the 
magnetic  field  lines  [figure  from  Keskinen  and  Ossakow,  1982]. 
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Fig.  8  —  Measured  electron  temperature  profile  (T{J\  solid  line)  and  corresponding 
electron  temperature  profile  obtained  from  the  quasilinear  theory  (T^,  dashed  line) 
for  three  different  electric  fields  (85)  mV/m,  upper  panel;  65  mV/m,  middle  panel; 
and  40  mV/m,  lower  panel).  The  CRA  neutral  temperature  Tn  has  also  been  plotted 
for  comparison  [from  St. -Maurice  et  at.,  1981 J. 
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MCDONNEL  DOUCLAS  CORPORATION 
5301  BOLSA  AVENUE 
HUNTINGTON  BEACH,  CA  92647 
01 CY  ATTN  N.  HARRIS 
OICY  ATTN  J.  MOULE 
01 CY  ATTN  GEORGE  MROZ 
01  CY  ATTN  W.  OLSON 
01  CY  ATTN  R.W.  HALPRIN 
01 CY  ATTN  TECHNICAL  LIBRARY  SERVICES 


PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  3027 
BELLEVUE,  WA  98009 

01 CY  ATTN  E.J.  FREMOUW 

PHYSICAL  DYNAMICS,  INC. 
P.O.  BOX  10367 
OAKLAND,  CA  94610 
ATTN  A.  THOMSON 


MISSION  RESEARCH  CORPORATION 
735  STATE  STREET 
SANTA  BARBARA,  CA  93101 


OICY 

ATTN 

P.  FISCHER 

OICY 

ATTN 

W.F.  CREVIER 

OICY 

ATTN 

STEVEN  L.  GUTSCHE 

OICY 

ATTN 

D.  SAPPENFIELD 

OICY 

ATTN 

R.  BOGUSCH 

OICY 

ATTO 

R.  HENDRICK 

OICY 

ATTO 

RALPH  KILB 

OICY 

ATTO 

DAVE  SOWLE 

OICY 

ATTO 

F.  FAJEN 

OICY 

ATTO 

M.  SCHEIBE 

OICY 

ATTO 

CONRAD  L.  LONGMIRE 

OICY 

ATTO 

WARREN  A.  SCKLUETER 

MITRE  CORPORATION,  THE 
P.O.  BOX  208 
BEDFORD,  MA  01730 

01CY  ATIN  JOHN  MORGAN  STERN 
01 CY  ATTN  G.  HARDING 
01CY  ATTO  C.E.  CALLAHAN 


R  &  D  ASSOCIATES 

P.O.  BOX  9695 

MARINA  DEL  REY,  CA  90291 


OICY 

ATTO 

FORREST  GILMORE 

OICY 

ATTO 

BRYAN  GABBARD 

OICY 

ATTN 

WILLIAM  B.  WRIGHT, 

OICY 

ATTN 

ROBERT  F.  LELEVIER 

OICY 

ATTO 

WILLIAM  J.  KARZAS 

OICY 

ATTO 

H.  ORY 

OICY 

ATTO 

C.  MACDONALD 

OICY 

ATTN 

R.  TURCO 

RAND  CORPORATION,  THE 
1700  MAIN  STREET 
SANTA  MONICA,  CA  90406 
01  CY  ATTN  CULLEN  CRAIN 
01 CY  ATTN  ED  BEDROZIAN 

RAYTHEON  CO. 

528  BOSTON  POST  ROAD 
SUDBURY,  MA  01776 

01CY  ATTN  BARBARA  ADAMS 


MITRE  CORP 

WESTGATE  RESEARCH  PARK 
1820  DOLLY  MADISON  BLVD 
MCLEAN,  VA  22101 
01 CY  ATTN  W.  HALL 
01 CY  ATTN  W.  FOSTER 

PACIFIC-SIERRA  RESEARCH  CORP 
1456  CLOVERFIELD  BLVD. 

SANTA  MONICA,  CA  90404 
OICY  ATTN  E.C.  FIELD,  JR. 

PENNSYLVANIA  STATE  UNIVERSITY 
IONOSPHERE  RESEARCH  LAB 
318  ELECTRICAL  ENGINEERING  EAST 
UNIVERSITY  PARK,  PA  16802 
(HO  CLASS  TO  THIS  ADDRESS) 

OICY  ATTN  IONOSPHERIC  RESEARCH  LAB 

PHOTOMETRICS,  INC. 

442  MARRETT  ROAD 
LEXINGTON,  MA  02173 

OICY  ATTN  IRVING  L.  KOFSKY 


RIVERSIDE  RESEARCH  INSTITUTE 
80  WEST  END  AVENUE 
NEW  YORK,  NY  10023 

OICY  ATTN  VINCE  TRAPANI 

SCIENCE  APPLICATIONS,  INC. 
P.O.  BOX  2351 
LA  JOLLA,  CA  92038 

OICY  ATTN  LEWIS  M.  LINSON* 
OICY  ATTN  DANIEL  A.  HAMLIN 
OICY  ATTN  E.  FRIEMAN 
OICY  ATTN  E.A.  STRAKER 
OICY  ATTN  CURTIS  A.  SMITH 
OICY  ATTN  JACK  MCDOUGALL 

SCIENCE  APPLICATIONS,  INC 
1710  COODRIDCE  DR. 

MCLEAN,  VA  22102 
ATTN:  J.  COCKAYNE 
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SRi  INTERNATIONAL 
333  RAVENSWOOD  AVENUE 
MENLO  PARK,  CA  94025 


01CY 

ATTN 

DONALD  NEILSON 

0  ICY 

ATTN 

ALAN  BURNS 

01CY 

ATTN 

G.  SMITH 

01CY 

ATTN 

L.L.  COBB 

01CY 

ATTN 

DAVID  A.  JOHNSON 

01CY 

ATTN 

WALTER  G.  CHESNUT 

01CY 

ATTN 

CHARLES  L.  RINO 

0  ICY 

ATTN 

WALTER  JAYE 

0  ICY 

ATTN 

M.  BARON 

01CY 

ATTN 

RAY  L.  LEAD AB RAND 

01CY 

ATTN 

G.  CARPENTER 

01CY 

ATTN 

G.  PRICE 

01CY 

ATTN 

J.  PETERSON 

0  ICY 

ATTN 

R.  HAKE,  JR. 

01CY 

ATTN 

V.  GONZALES 

01CY 

ATTN 

D.  MCDANIEL 

STEWART  RADIANCE  LABORATORY 
UTAH  STATE  UNIVERSITY 
1  DE  ANGELO  DRIVE 
BEDFORD,  MA  01730 
01 CY  ATTN  J.  ULWICK 

TECHNOLOGY  INTERNATIONAL  CORP 
75  WIGGINS  AVENUE 
BEDFORD,  HA  01730 

01CY  ATTN  W.P.  BOQUIST 

TRV.  DEFENSE  &  SPACE  SYS  GROUP 
ONE  SPACE  PARK 
REDONDO  BEACH,  CA  90278 
01CY  ATTN  R,  K.  PLEBUCH 
01CY  ATTN  S.  ALTSCHULER 
01CY  ATTN  D.  DEE 

VISIDYNE 

SOUTH  BEDFORD  STREET 
BURLINGTON,  MASS  01803 
01CY  ATTN  W.  REIDY 
0 ICY  ATTN  J.  CARPENTER 
0 ICY  ATTN  C.  HUMPHREY 
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IONOSPHERIC  MODELING  DISTRIBUTION  LIST 
(UNCLASSIFIED  ONLY) 


PLEASE  DISTRIBUTE  ONE  COPY  TO  EACH  OF  THE  FOLLOWING  PEOPLE: 


NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  D.C.  20375 
DR.  P.  MANGE  -  CODE  4101 
DR.  R.  MEIER  -  CODE  4141 
DR.  E.  SZUSZCZEWICZ  -  CODE  4187 
DR.  J.  GOODMAN  -  CODE  4180 
Dr.  R.  RODRIGUEZ  -  CODE  4187 

A.F.  GEOPHYSICS  LABORATORY 
L.G.  HANSCOM  FIELD 
BEDFORD,  HA  01730 
DR.  T.  ELKINS 
DR.  W.  SWIDER 
MRS.  R.  SAGALYN 
DR.  J.M.  FORBES 
DR.  T.J.  KENESHEA 
DR.  J.  AARONS 
DR.  H.  CARLSON 
DR.  J.  JASPERSE 

CORNELL  UNIVERSITY 
ITHACA,  NY  14850 
DR.  W.E.  SWARTZ 
DR.  R.  SUDAN 
DR.  D.  FARLEY 
DR.  M.  KELLEY 

HARVARD  UNIVERSITY 
HARVARD  SQUARE 
CAMBRIDGE,  MA  02138 
DR.  M.B.  McELROY 
DR.  R.  LINDZEN 

INSTITUTE  FOR  DEFENSE  ANALYSIS 
400  ARMY/ NAVY  DRIVE 
ARLINGTON,  VA  22202 
DR.  E.  BAUER 

MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
PLASMA  FUSION  CENTER 
LIBRARY,  NW16-262 
CAMBRIDGE,  MA  02139 


RASA 

GODDARD  SPACE  FLIGHT  CENTER 
GREENBELT,  MD  20771 
DR.  S.  CHANDRA 
DR.  K.  MAEDA 
DR.  R.F.  BENSON 

NATIONAL  TECHNICAL  INFORMATION  CENTER 
CAMERON  STATION 
ALEXANDRIA,  VA  22314 
12CY  ATTN  TC 

COMMANDER 

NAVAL  AIR  SYSTEMS  COMMAND 
DEPARTMENT  OF  THE  NAVY 
WASHINGTON,  D.C.  20360 
DR.  T.  CZUBA 

COMMANDER 

NAVAL  OCEAN  SYSTEMS  CENTER 
SAN  DIEGO.  CA  92152 

MR.  R.  ROSE  -  CODE  5321 

NOAA 

DIRECTOR  OF  SPACE  AND  ENVIRONMENTAL 
LABORATORY 
BOULDER,  CO  80302 
DR.  A.  GLENN  JEAN 
DR.  G.W.  ADAMS 
DR.  D.N.  ANDERSON 
DR.  K.  DAVIES 
DR.  R.  F.  DONNELLY 

OFFICE  OF  NAVAL  RESEARCH 
800  NORTH  QUINCY  STREET 
ARLINGTON,  VA  22217 
DR.  G.  JOINER 

PENNSYLVANIA  STATE  UNIVERSITY 
UNIVERSITY  PARK,  PA  16802 
DR.  J.S.  NISBET 
DR.  P.R.  ROHRBAUGH 
DR.  L.A.  CARPENTER 
DR.  M.  LEE 
DR.  R.  DIVANY 
DR.  P.  BENNETT 
DR.  F.  KLEVANS 


PRINCETON  UNIVERSITY 
PLASMA  PHYSICS  LABORATORY 
PRINCETON,  NJ  08540 
DR.  F.  PERKINS 

SCIENCE  APPLICATIONS,  INC. 
1150  PROSPECT  PLAZA 
LA  JOLLA,  CA  92037 
DR.  D.A.  HAMLIN 
DR.  L.  LINSON 
DR.  E.  FRIEMAN 

STANFORD  UNIVERSITY 
STANFORD,  CA  94305 
DR.  P.M.  BANKS 


UNIVERSITY  OF  PITTSBURGH 
PITTSBURGH,  PA  15213 
DR.  N.  ZABUSKY 
DR.  M.  BIOND1 
DR.  E.  OVERMAN 

UTAH  STATE  UNIVERSITY 
4TH  AND  8TH  STREETS 
LOGAN,  UTAH  64322 
DR.  R.  HARRIS 
DR.  K.  BAKER 
DR.  R.  SCHUNK 


U.S.  ARMY  ABERDEEN  RESEARCH 
AND  DEVELOPMENT  CENTER 
BALLISTIC  RESEARCH  LABORATORY 
ABERDEEN,  MD 

DR.  J.  HEIMERL 


UNIVERSITY  OF  CALIFORNIA, 
BERKELEY 

BERKELEY,  CA  94720 
DR.  M.  HUDSON 


UNIVERSITY  OF  CALIFORNIA 
LOS  ALAMOS  SCIENTIFIC  LABORATORY 
J-10,  MS-664 
LOS  ALAMOS,  NM  87545 
M.  PONGRATZ 
D.  SIMONS 
G.  BARASCH 
L.  DUNCAN 
P.  BERNHARDT 


UNIVERSITY  OF  CALIFORNIA, 
LOS  ANGELES 
405  HILLGARD  AVENUE 
LOS  ANGELES,  CA  90024 
DR.  F.V.  CORONITI 
DR.  C.  KENNEL 
DR.  A.Y.  WONG 


UNIVERSITY  OF  MARYLAND 
COLLEGE  PARK,  MD  20740 
DR.  K.  PAPADOPOULOS 
DR.  E.  OTT 
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